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The paper presents Monte-Carlo simulated lateral dis- 
tribution functions for electrons of EPC developing in lead, 
at superhigh energies (.1-1 PeV) for depths t^60 c.u. 

A t=1 t Q c.u. The higher moment characteristics, i.e. varia- 
tion, asymmetry, excess, are presented along with analyti- 
cal solutions for the same characteristics at fixed obser- 
vation level calculated to theory approximations A and B 
by using numerical inversion of the Laplace transformation. 
The conclusion is made of a complex, usually non-Gaussian 
shape of the function of the particle number distribution 
within a circle of given radius at fixed depth. 

To analyse experimental data obtained by an X-ray 
emulsion chamber technique, the detailed information on 
mean EPC characteristics and their fluctuations in dense 
media is necessary. Therefore, in the high energy region 
The Landau-Pomeranchuk effect should be regarded, in the 
low energy region ionization losses and the Compton effect 
should be allowed for, and it is necessary to describe ra- 
ther correctly the scattering processes. The problem being 
extremely complex requires the Monte-Carlo method to be 
used. For a simplified problem(one-dimensional theory ap- 
proximations A and B, i.e. consideration of the number of 
particles with energy higher than the given energy), it is 
possible now to analyse analytically the problem of the 
distribution function for the particle number at fixed 
depth using higher moments, unlike t.o /I/, where only the 
second moment has been used. This approach allows qualita- 
tive comparison with the Monte-Carlo results and is of 
interest in itself as^an approbation of the technique (e.g. 
when using the latter to analyse quark-gluon cascades). 

1 . A cascade from a primary of energy E 0 was Monte- 
Carlo simulated up to energies about E=10~2 -io~3e 0 , then 
each branch of 130-150 branches of the cascade was ended 
by statistically with a mean cascade from a previously 
calculated data bank. In calculations the above mentioned 
processes were regarded, and scattering at each segment 
of the free path was considered in the Fermi approxima- 
tion/2/. An estimate by the method of standards assuming 
the initial distribution to be non-Gaussian/34, showed 
that, in the range (1/2v2)t max , where tmax is the depth 
of a mean shower maximumCtmax^Oto) , a statistics of 100 
events yielded an error (standard) of the mean to be 2-3%, 
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of variation *■'10%, of asymmetry 20-30$, and of excess 
40-80$. Calculations were made on a net of 14 radii Cl -400 
um) given uniformly in logarithm of spatial variable, (lgr), 
that made it feasible to plot the differential P(E 0 ,r,t) 
lateral distribution function: the particle number sdb dis- 
tance r from the axis at depth t. The obtained differential 
functions were integrated and approximated. 

Pigs.1 and 2 compares our calculations (see the solid 
lines) with the experiment/4/ for the differential and in- 
tegral distribution functions. Thw experimental restilts are 
presented for front-side film of a type C-chambers woth 
delution factor d=1.18 leading to a decrease of the particle 
number in cascade with increasing depth. 

Pig. 3 shows the results for 8 values of primary Y~ . 

The bundles of curves are parametrized by the product E 0 R» 
The point in Pig. 3 are for the results/5/. *his figure il- 
lustrates violation of the core approximation of the Lan da u— 
Pomeranchuk effect. The investigation of fluctuations of the 
particle number within a circle of given R substantiates the 
general picture of evolution of the density distribution 
function/6/ - at t^1/2 t nw-r and t>3/2 tmax» the distribu- 
tions are narrow-peaked and skewed towards the small par- 
ticle number (to the left), at 1/2 t moy t ^4 3/ 2 tmax these 
distributions shift to the right transforming through sym- 
metric distributions into asymmetric ones with a flat maxi- 
mum(or several maxima) , Gaussian distributions in the region 
of tmax are rather an exclusion than a rule. The results 


tinder consideration give a considerable qualitative clarifi- 
cation of the earlier data and allow investigation of the 
dependence on R or, what is the same, on the efficient 


um(b). The influence of the Landau-Pomeranchuk effect can 
be clearly seen in Pig. 4 - the curve minima 0 !” E 0 10 TeV 
(see Pig. 4a); degradiag(see Pig. 4b) with increasing R cor- 
responds to a decrease of efficient threshold energy, i.e. 
a weaker influence of the landay-Pomeranchuk effect on a 
cascade. 


2. To analyse analytically the set of equations for 
higher moments of the function of distribution over the par- 
ticle number in approximation A, it is written in the simi- 
lar form as in /I/. The solution of this set of equations 
can be obtained to approximation A for homogeneous cross 
sections and any moment using the Mellin transformation in 
energy E 0 /l/. Por the first, second, third etc. moments of 
the functions of distribution over the particle number in a 
shower Induced by a primary electron the solutions can be 
represented as follows /7/: 

~T~j / _ ,\ j ,fls { 0 +A T 


f(c ^0 ~ J AA (^/e) * 

- t ^ M 

Hi (tTn/kv) "si?! 
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In the cascade theory approximation B similar formulae 
for higher moments can be obtained. These integrals for a 
moderate values C ti =1 ,2,3,4) can be calculated using the 
numerical Laplace inversion. These caleulational results are 
listed in Table 1 and show qualitative agreement with the 
simulation results. Around the shower maximum we can per- 
form the detailed analysis of higher moments and structure 
of the distribution function at E 0 /E(E 0 /i&) . Using the 

saddle point method (S -1) it can be shown that in the 
region of shower maxim&fr * 

Van ^ A 2 //<£ 7 &/£ 7 « — —7 

i.e. Var decreased slowly, as /«•&/£ , with increasing ratio 
E 0 /E/8/. Asymmet ry and excess increases _ ■ 

A x ^ Aj f&aTF , Ex ~ A 4 % 

This behaviour of higher moments contradicts the assump- 
tion of the Gaussian shape of the distribution functions for 
the particle number. 

The investigation showed that fluctuations of the dist- 
ribution functions for the particle number Within a 
circle of fixed radius and those of the function of the num- 
ber of particles with energy higher than the given energy 
are of similar character, oscillations of higher moments 
show evidence of a complex structure of the distribution 
functions and indicate that Gaussian distribution is rather 
an exclusion than a rule. 
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Table 1 


t 

U Var 2 A 

f 

Var 2 

A s 

9 

T 

Var 

A s 

2 

64.8 .25 1.0 

30.1 

.17 

•1.92 

.12+2 

*25 

.99 

4 

415 .14 .3 

110.2 

-.It 

-.12 

.13+4 

.15 

.52 

6 

982 .7-1 -.45 

157.9 

.48-1 

-1.1 

.45+4 

.94-1 

-.4-2 

8 

1333.31-1-.78 

145 

.33-1 

-1.3 

— — 

— 

-.2+1 

10 

1272.17-1-1.26 

91.9 

.76-1 

.76-1 

.12+5 

.24-1 

-3.34 
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Table 1 (continued) 

t 5 Var 2 A„ 5 Var 2 A a I Var 2 17 

8 p p 

14 600 . 9-1 .6 24.5 .49 1.75 . 95+4 . 15-1 .87 
18 167 . 43-1 .91 4.44 2.29 4.99 . 40+4 .10 1.77 
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